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We consider a system consisting of a 2D network of links between Majorana fermions on super- 
conducting islands. We show that the fermionic Hamiltonian modeling this system gives rise to 
Kitaev's toric code in fourth-order perturbation theory. By using a Jordan- Wigner transformation 
we can map the model onto a family of signed 2D Ising models in a transverse field where the 
signs, ferromagnetic (FM) or anti-ferromagnetic (AFM), are determined by additional gauge bits. 
Our mapping allows an understanding of the non-perturbative regime and the phase transition to a 
non-topological phase. We discuss the physics behind a possible implementation of this model and 
argue how it can be used for topological quantum computation. 

PACS numbers: 03.67.-a, 03.65.Vf, 03.67.Lx 



Kitaev's well-known toric code 0, is a toy model 
Hamiltonian which demonstrates the concept of topolog- 
ical order in two dimensions. One can imagine storing 
a qubit in this ground space of this model such that, at 
low temperature T compared to the gap, dephasing of 
such qubit is exponentially suppressed with growing lat- 
tice size. More general constructions allow the encoding 
of many qubits in the ground space and the topological 
implementation of a cNOT and Hadamard gate by means 
of Hamiltonian or code deformation 




FIG. 1: (Color online) (a) Fermionic model studied in this 
paper. Each island (light gray square) has four Majorana 
fermions (yellow dots) labeled as a, b, c, d. 50% occupancy is 
favored for these two fermionic modes, as expressed by the 
parity constraint in Eq. fTJ. A weaker quadratic interaction 
exists between Majorana fermions on islands i and j along 
diagonal links Vi =fJ ,±zj =tJ ,±x, Eq. ([2]). (b) Zoomed view of 
a single island. The Majorana wire has a C shape (black 
line) in order to be able to tune the overlap between the 
a- and c-Majorana and thus implement an X-gate. A Z- 
gate/measurement is implemented by increasing the ratio of 
E' c /E'j, see main text. 



Majorana bound states at their ends IB" IB], we believe 
that our model may provide a viable route to the re- 
alization of topological quantum computation. The in- 
terest in Majorana fermion wires is partially motivated 
by their fermionic-parity protected ground space degen- 
eracy which allows parity protected quantum computa- 
tion [nHIl and braiding in networks of nanowires [20l — 
[23}. The advantage of the approach advocated in this 
paper is that the protection is fully topological and no 
longer based on fermionic-parity conservation. The idea 
of engineering a topologically-ordered Hamiltonian using 
Josephson-junction arrays has been explored mostly in 
the work of Ioffe et al, see e.g. [2fjj . 

We consider the following fermionic Hamiltonian H = 

H + V where H = Y1^ L ^0 an d * labels the square 
islands in Fig. Q] The lattice in Fig. [IJa) has periodic 
boundary conditions in both directions (see the Appendix 
for a discussion of the model with open boundaries). 
Each Hq acts on two fermionic modes or four Majorana 
modes as 



Tji \ i i i i 

H = - Ac a c bC c C d - 



(1) 



Further, we have V — A J2 i< j Vi.j (with A > 0) where 
i,j represents the interaction between two Majorana 
fermions on adjacent islands i and j, i.e., Vi =tl ±z j =J u±x 
for a plaquette [i equals 

V - - ±ic^ ± c^ +i V - - =+icl Jt ~ i ( Ji ~ i 

'IM+z,fi—x J - M 'b °c > v iJ,—z,fi—x — -*- lu a u d ' 



(l — Z,fJ,-\-X 



±ict +£ c^, V, +w = ±ic% +S c%+*. (2) 



In this paper, we investigate how one could arrive at 
the toric code Hamiltonian starting from a realistic model 
of interacting fermions. Kitaev has shown how the toric 
code emerges in fourth-order perturbation theory from 
the so-called honeycomb model [|[ (see also [§]). Given 
the recent interest in making proximity-coupled semi- 
conducting nanowires which support weakly-interacting 



All link operators V^±z^±x mutually commute. The 
± signs of these terms will be fixed according to the con- 
sistent orientation of the plaquettes in Fig. [TJ i.e., the 
link on the top-left of a white plaquette /i represents the 
interaction Vp+s^s = ic%~ x c% +z . Physically, the signs 
depend on microscopic detail and thus we assume them 
to be random but fixed. We can find an extensive set of 
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FIG. 2: (Color online) (a) Toric code on a L x L lattice with 
2L 2 qubits on vertices. The Hamiltonian is a sum over white 



and grey plaquette operators 



Zi \i-\-z Z \i — z ^(i-fi^-i ■ 



(b) 



When the hatched plaquettes are omitted in the toric code 
Hamiltonian, the ground space degeneracy increases from 4 to 
8. The logical operators of this additional qubit are depicted 
by the blue and red lines. We can call this a white hole qubit 
as it is obtained by making a hole in the lattice which is 
centered around a white plaquette. 



operators which commute with all terms of H and which 
all mutually commute. These are, first of all, weight- 
8 fermionic plaquette operators {C/J,C™} (where g (w) 
stands for gray (white) plaquettes fx) which are the prod- 
uct of four link operators around a plaquette: 



°d L a L c L 6 L a 



-zji-xu-x u+z 
c d c b L c 



(3) 



Secondly, the torus has two homologically non-trivial 
closed loops 71,72 and the loop operators C 7l/2 = 
II(j ,j)e 7l/2 Vij commute with all link Vij and island op- 
erators Hq. Thus the Hamiltonian is block-diagonal 
with respect to subspaces ('sectors') characterized by the 
eigenvalues C 9 J W — ±1 and C 7l/2 = ±1. 

We analyze the model in the perturbative regime where 
A ^> A, see e.g. [§]. The ground space of H at A = is 
characterized by {c a c\d c c d = +1}, and thus the ground 
space on each island is a two-dimensional subspace, a 
qubit. One can define the logical X and Z operator on 
this island qubit as 



Xi 



(4) 



Let P_ = 2 2L2 ltftl(I + 4444) be the projector onto 
this 2 2L -dimensional unperturbed ground space. The 
first non-trivial term in the perturbative expansion oc- 
curs in 4 th order (see the Appendix), i.e. 
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where A.. 



i .e., the plaque- 



tte terms of the toric code in Fig. [2] [28). Note that 
P-Cf/ W P- = A 9 J W and hence the four-dimensional toric 



code ground space of H when A ^> A lies in the {C, 
+1} sector. 
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FIG. 3: (Color online) The gauge bits a set the Ising inter- 
actions to FM (black edges) except for an AFM (red edges) 
loop around the torus. This AFM boundary will be felt in the 
ferromagnetic (FM) phase, but not in the paramagnetic (PM) 
phase of the model, leading to degeneracy. A loop operator 
C~, in the fermionic model becomes a product of Ising edges 
which winds around the torus. 



Let us consider how the topological phase extends to the 
regime where f = 0(1). We can get insight for this 
regime by considering higher-order terms in the pertur- 
bative expansion (see e.g. [29} for such expansion for 
the honeycomb model), see the Appendix. This calcula- 
tion shows that the ground space degeneracy of the toric 
code is broken only in (2L) th -order perturbation theory 
and we expect that the topological phase will destabilize 
via this mechanism. 

To consolidate this picture we map our model via a 
Jordan- Wigner (JW) transformation, onto a family of 
signed transverse field (TF) Ising models on a 2D square 
lattice, see Appendix. The JW mapping preserves the 
locality of the interactions, and all reductions are ex- 
plicit and rigorous. Previous work 27[ has shown how to 
map an Ising gauge theory in a transverse field (assuming 
open boundary conditions) by a Jordan- Wigner transfor- 
mation onto the Majorana fermion model discussed here. 
Our mapping clarifies the nature of the topological phase 
transition and the parameter values for which it will oc- 
cur. In this analysis we restrict ourselves to the sector 
{C™ = +1}. Note that when A = the ground state of 
H is unique since we are in a state of fixed link parity 
{V hJ = -1}. This state lies in the {C a J w = +1} sec- 
tor since each CfJ w is a product of four link operators. 
As we saw above, at A -C A, the ground space also lies 
in the {C g J w = +1} sector. The mapping decouples our 
fermionic model into a set of transverse field Ising models 
H(a) = -XJ2i i <?i,jSiS- - A^ Sf with the condition 



^-(i.j)e white = 1 so that the sign of the Ising inter- 
actions is determined by gauge bits a^j € {—1,1} as- 
sociated with the edges. The gauge condition expresses 
the fact that the white plaquettes are never frustrated, 
i.e. = +1, but a gray plaquette is frustrated when 

c?, = n, 



(ij)egray IJ. a i,j 



— 1. The spectrum of H{a) 
solely depends on the frustration of the Ising interac- 
tions and the presence of domain walls or homologically 
non-trivial loops, see Fig. [3] We anticipate the follow- 
ing spectrum, see Fig. @] At both ends of the parameter 
region (A or A = 0), the groundspace lies in the un- 
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FIG. 4: Sketch of the spectrum of system as a function of 
A/A. For small A, the system is in a topological state with 
a four-fold ground state degeneracy on the torus. The first 
excited states for small A are Ising models with frustration 
as determined by the gauge bits. All these models are de- 
generate for A = 0; and the degeneracy lifts in forth order 
perturbation theory in A, see Eq. ((5j). The phase transition 
to a state without topological order happens at the transition 
point (A/A) c of the unfmstrated Ising model. At this transi- 
tion the gap of the Ising model closes and the degeneracy of 
the topological states vanishes. 



frustrated TF Ising model sector. We have numerically 
confirmed this for small lattice sizes for the entire param- 
eter regime. Because of the symmetry between gray and 
white plaquettes, this finding also motivates the choice 
for {C™ = +1} as the ground sector. Fig.[3]depicts a con- 
figuration a which represents a homologically non-trivial 
loop; all plaquettes are unfrustrated, but an Ising model 
with such an AFM sign pattern will contain a domain 
wall of length at least L where bonds are not satisfied. 
The topological phase is identified as the paramagnetic 
(PM) phase (S z ) = in the TF Ising models. In this 
phase the ground spaces of the Hamiltonians H(a) with 
unfrustrated configurations a with or without the 2 non- 
trivial AFM loops are approximately degenerate: this is 
the topological degeneracy whose splitting we expect to 
scale as exp(— where £ is the correlation length of 
the TF ferromagnetic Ising model. We expect the ef- 
fective gap A e ff above the degenerate ground-space to 
monotonically increase before we reach the second-order 
phase transition of the TF ferromagnetic Isin g m odel, 
which is known to occur around (■£■) ~ 0-33 [30(. El- 
ementary excitations of the toric code with = +1 
for two gray plaquettes correspond to ground-states 
of TF Ising models with frustration at those particular 
gray squares. 



We can consider the effect of additional perturbations. 
One would expect quadratic Majorana fermion pertur- 
bations of strength e on each island due the coupling of 
the wire ends on an island. If a qubit is encoded in a pair 
of Majorana wires, such coupling leads to an energy-level 



splitting of the qubit state of strength e. Here the advan- 
tage of our topological encoding becomes clear. Consider 
a perturbation U oc eic^ ~ x c%~ x where /i is some white 
plaquette (by symmetry other perturbations would act 
similarly) and e -C A. On the groundspace of H , the 
perturbation acts as a local term oc eX M _^, hence we ex- 
pect that the topological degeneracy is preserved up to 
some critical strength (3^)0 where A e g is the gap above 
the degenerate ground-space in the topological phase. In 
practice, we expect these quadratic perturbations to be 
random (and weak) and hence they could be beneficial 
in stabilizing the topological quantum memory at finite 
temperature by limiting the diffusion of anyons (see e.g. 

We analyze the possibility of implementing the model 
presented above in superconducting-semiconducting het- 
erostructures. Putting a semiconducting nanowire such 
as InAs or InSb with strong spin-orbit interaction on top 
of a conventional superconductor subject to a sufficiently 
strong magnetic field leads to a pair of Majorana modes 
located at the ends of the nanowire HlQ. We focus 
on a 2D array of superconducting islands each supporting 
two nanowires leading to four unpaired Majorana modes, 
see Fig. [TJ The product of the four Majorana operators 
V 1 = c l a c\c l c c l d is fixed by the parity V % = (— lV" of the 



number of electrons rij on the i-th island [24|, |3l| . Two 



Majorana modes c x and c y on different superconducting 
islands i and j interact with each other via the anomalous 
Josephson interaction 



jV itj cos[(c/>j 



<M/2], 



(6) 



where fa denotes the superconducting phase on the i-th 
island and Tij is proportional to the probability ampli- 
tude for tunneling a single electron/hole across the link 
i, j from mode c l x to c? y between the islands i and j 
Along the lines of Ref. , we shunt each superconduct- 
ing island with a strong Josephson coupling Ej to a com- 
mon ground superconductor. This Josephson coupling 
fixes the superconducting phases fa to a common value 
due to large fluctuations of the number of electrons (in 
units of two) on and off the island. Note that this way 
of freezing out the bosonic degrees of freedom due to the 
superconductivity is different from the one discussed in 
Ref. 2_4| which proposes a large charging energy which 
fixes the number of electrons by delocalizing the super- 
conducting phase completely. Even though the charge is 
strongly fluctuating, the fermion parity V 1 remains con- 
served. In the ideal case when all the Josephson couplings 
are equal r^- = A, the anomalous Josephson interaction 
Hj implements V. Of course, V in general will not have 
the orientation indicated in Fig. [1] but if we work with 
a lattice with open boundary conditions (corresponding 
to the surface code [32|]), any sign pattern of the link 
interactions simply picks out a topological sector with 
a corresponding pattern of eigenvalues A^ = ±1 as the 
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ground-space [331 ] . see our discussion in the Appendix. 
From the coding perspective it is well known that topo- 
logical computation can proceed not just in the trivial 
syndrome (all eigenvalues of = +1) sector but also in 
some non-trivial syndrome sector. 

Next, we discuss the effect of a capacitive coupling 



H C = E c y^X n i - n ind) 2 



(7) 



to the ground plate with the capacitive energy Ec — 
e 2 /2C where C is the capacitance of the island with 
respect to ground and ni n( j a constant which is due to 
gate voltages. For simplicity of notation, we have as- 
sumed all the capacitances C and offset charges nj n( j to 
be equal. In the regime with A < Eq <C Ej, the capac- 
itive coupling introduces phase-slips through the strong 
Josephson junctions and thus an energy difference be- 
tween states with different fermion parities. This leads 

to Eq. © with A oc Ec /4 Ej /4 cos(Ttn ind )e-^ SEj/Ec 
The sign of A depends on the value of ni n d can thus 
be tuned in principle. A residual interaction between 
two Majorana modes c l x and c l y with strength e is due 
to overlap of the wave functions of the Majorana bound 
states. However, because the states are localized, this 
coupling can be made exponentially small by keeping the 
modes sufficiently far apart from each other. Having suf- 
ficiently strong tunneling coupling between the Majorana 
fermions along the links, a value A ~ 200 mK seem real- 
istic as the bare gap of the Majorana wire is likely of the 
order of a few K [20|. Choosing the Josephson energy 
to be Ej ~ 10 K and a capacitive energy Ec — 5K, we 
obtain Ej A > A as required . The resulting optimal 
value of A e g will be of the order of A ~ 200 mK. 

It is possible to tune the tunnel coupling T^j by chang- 
ing the tunneling barrier between island i and j by a 
nearby gate. This is a mechanism through which we can 
create holes in the lattice. In practice one can work with 
a lattice of superconducting islands which represents the 
surface code with open boundary conditions, encoding 
one logical qubit. One can apply the ideas of the surface 
code architecture if (i) one can make (and move) gray 
and white holes of arbitrary size in this surface and (ii) 
one can locally measure and Zi and prepare JQ and 
Zi eigenstates. The preparation of certain 1-qubit an- 
cillas can then be used to achieve universality [35|. An 
example of a white hole qubit is depicted in Fig. [2h). 
Such a white hole can be obtained by cutting off the four 
black links surrrounding the center white plaquette such 
that the hatched plaquette terms no longer appear in the 
effective Hamiltonian. Moving such a white hole could be 
done by adiabatically changing the strength of Majorana 
links in order to turn links on and off. The operations (ii) 
can be implemented using the setup of Fig.QJb). Instead 
of a single superconducting island, each site consist in 
fact of two islands with two Majorana modes each. Most 
of the time, these islands are coupled to each other via 



a strong Josephson coupling E'j E' c ,Ec (see Fig. [T]) 
such that they essentially act as a single island such that 
all the discussion above applies unchanged. Increasing 
the ratio E' c /E'j turns on a magnetic field along the z- 
axis which can be used to implement rotations around 
this axis. Additionally, the measurement of Zi can be 
implemented by coupling one of the superconducting is- 
lands to a fermion parity meter [13] ■ Single qubit univer- 
sality is achieved by increasing the overlap of the a and c 
Majorana modes — by decreasing the length of the topo- 
logical trivial part of the Majorana wire indicated by the 
dashed — and thus effectively implementing a magnetic 
field along the x-axis. 

We acknowledge fruitful discussions with A. 
Akhmerov. DDV and FH are grateful for support 
from the Alexander von Humboldt foundation. 



Appendix 

We start by detailing some of the steps in the perturba- 
tive analysis. We defined P_ = 2~ 2L2 nff 1 (7 + 4444) 
as the projector onto the 2 2L -dimensional groundspace 
and P- + P+ = I. Let V± T = P±VP T and let G+ = 
P + Hq 1 P + where we have redefined Ho as Hq + 2AL 2 
such that its lowest-eigenvalue is 0. In the self-energy 
expansion all terms with odd number of perturbations V 

vanish. The second-order term contributes 

a term proportional to I whereas the fourth-order term 
equals 



V-+G+V++G+V++G+V+- = - 
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where A^ — 



Zi_ t +zX n 



u+zZn-iXf.-*, i.e., the plaquette 
terms of the toric code in Fig. [2] [28| . Note the difference 
with Kitaev's honeycomb model where the prefactor of 
these plaquette terms is instead of j^; this is be- 
cause the links in our model commute whereas the links 
in Kitaev's model mutually anti-commute on their com- 
mon qubit. All higher-order terms in the perturbative 
expansion consist of loops of links on the lattice and can 
be represented as products of plaquette operators A^, 
i.e. the terms are of the form — . . . A"i ■ ■ • A"* > 
all with the same negative sign, further stabilizing the 
ground space with {A^ = +1}. 

In (2L) th -order one obtains also terms proportional to 
A 1 = PqC~ ( Pq oc Tli ei (Ti where i is a product over islands 



through which the loop 7 goes and oi = Xi, Y% or Zi as 
defined in Eq. ((4]) depending on the loop 7 taking one of 
the three directions across island i. S ince A^ 1 and A.y 2 
commute with all {A^} and with each other, they must 
leave the ground space of the toric code invariant and be 
mutually commuting products of the 4 logical operators 
of the toric code qubits, i.e. they can be represented as 
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X1X2 and Z1Z2 where (Xi,Zi), i = 1,2 are the non- 
local logical operators of the 2 toric code qubits. Hence 
the presence of these terms in the effective Hamiltonian 
breaks the topological degeneracy. 

The overall mapping via the JW transformation goes 
as follows. We will denote the eigenvalues of the oper- 
ators CT, Eq. ([3]), as c™. Note, by the way, that there 
is one linear dependency between the plaquette opera- 
tors i.e. C a n = n M C^ = II^C™ where C a u is the prod- 
uct of all 8L 2 Majorana operators. First, we map our 
model via a JW transformation onto a model which for 
fixed eigenvalues {c™} consists of XX (strength A) and 
ZZ links (strength A), a square-octagon model on the 
left in Fig. |6l We consider this square-octagon model 
for {c™ = +1}, but we extend the state space to any 
state with ILjC™ = +1 [38|. In this extended state 
space, we use the additional symmetry of the square- 
octagon model to lay out a basis of Bell states on the 
islands. In the Bell basis, the square-octagon model re- 
duces to a set of transverse field Ising models where the 
sign of the Ising interactions is determined by an ad- 
ditional sign qubit. In principle this degree of freedom 
is present at every island; the sign qubits are the black 
and red dots in Fig. [6] However, these models are uni- 
tarily equivalent to ones in which associate a gauge bit 
(Ti.j e { — 1,1} with every edge and fix the gauge 

B(i j)G„hito n a i,3 — 1- Thus we obtain the transverse field 
Ising models are H(a) = AV , a,.,S;Sj - A^Sf 
where Hnj\ & whito ^Oi j — 1. Let us now explicitly show 
these steps. 

For the JW transformation we choose an order for 
the 8L 2 Majorana fermions ci,...c S l 2 - We will order 
the Majorana fermions around every white plaquette as 
in Fig. O i.e., we start with an arbitrary white pla- 
quette, choose this ordering of labels and go onto any 
next white plaquette and continue until we are done la- 
beling all Majorana fermions. Given this ordering the 
JW transformation maps T(c2j-i) = Z\ . . . Zi-\Xi and 
Y(c2i) = Z\...Zi—\Yi where the Pauli's here are of 
course unrelated to the Pauli's defined in Eq. ((4]). Under 
this JW transformation the links Vi = ^±z j =fJ ,±x around a 
white plaquette \i become 

— z,^i+x ^ X^ — zX^-\-x^ Vp,-\-z,p,-\-x ^ xZ^i— zY[i-\-x- 

The island operator of an island i which is north or south 
of a white plaquette becomes Hq —> — AZ^iZ^ where 
i, 1 and i, 2 label the two qubits on the island, see the 
blue links in Fig. [51 The island operator of an island 
i which is east or west of a white plaquette becomes 
Hq — > +AZi^Zi^- In this qubit representation the white 
plaquette operator equals — > —H^Q^Zi where \i now 
represents a square in Fig. [6] and i £ d/j, is the product 
over qubits at the corners of the square. Let the subspace 
with fixed eigenvalues {c™} be denoted as "%({c™}). We 
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FIG. 5: White plaquette where the anti-clockwise direction 
fixes the sign of the interaction V depicted by black links. 
For such a white plaquette we choose the Majorana fermion 
ordering indicated by the numbers 1,2. ..8 to do the Jordan- 
Wigner transformation. 

can represent the action of the weight-4 link operator on 
%{{c™}) as 

Vfl-\-X,[l+Z Y/l • zZ [i- x^fj. Z^'/X • X X fj,-\-zX f^ + x • (9) 

Note that on the space H({c^'}) the gray plaquette op- 
erator equals = c^_ £ _ £ Hi^g ll Xi. This implies that 
on each subspace "H({c™}) the action of the Hamiltonian 
can be best represented using a square-octagon lattice in 
which qubits live on the vertices, see Fig. [5] The white 
plaquettes have been transformed into squares and the 
gray plaquettes into octagons. The islands have been 
stretched in the horizontal and vertical directions and 
are represented by horizontal and vertical links equal 
to ±AZ itl Zi.2 (—1 for vertical links, +1 for horizon- 
tal links). On this lattice, diagonal links between ver- 
tices i and j are all XXiXj except the diagonal links 
Xc^X^+zX^+x for the white squares labeled by /i. In 
the remaining analysis we choose {c™ = +1}, but all 
further steps can be done for other {c™} [39j |. 



In order to use the square-octagon model, we extend 
H({c™ = +1}) into the (almost) full state space, obeying 
only the parity condition n^C™ = +1. At the end of this 
section, we show why this procedure only adds degenera- 
cies to the Hamiltonian, not altering the eigenspectrum. 
This parity constraint n M C™ = C a n = +1 is, after the 
JW transformation, equivalent to (— l) L YliZ; L = +1. A 
new symmetry then becomes apparent, namely the in- 
teractions of the square-octagon model commute with 
vertical and horizontal link operators Kj = Xj^Xj^- 
Note that by an even number of local unitary X rota- 
tions we can rotate all horizontal and vertical links to be 
— AZj \Zj 2- From now on we take L even for simplic- 
ity so that the parity constraint implies that YliZ L = 1. 
A basis for the state space can then be constructed by 
laying out Bell states \^> st ) = |s,i) on the horizontal 
and vertical links between qubits. Here s represents the 



6 



sign qubit or eigenvalue of XX, i.e., \s) corresponds to 
_^(| 00 ) + (-l) s |ll)) and ^(|01) + (-1)110)). Qubit t 
represents the eigenvalue of ZZ, i.e., \t) corresponds to 
(Ix <g) X|)^=(|00) ± 1 11)). This choice of basis and local 
unitary rotation to sign qubits s and Ising qubits t allows 
one to write the Hamiltonian on the square-octagon lat- 
tice as a transverse field Ising model on the t qubits on 
a square lattice where the additional sign qubits live at 
every site, see Fig. [51 

In order to obtain this model we lay out the Bell states 
such that any XX link in the model acts between qubit 1 
of one horizontal (resp. vertical) Bell state and qubit 2 of 
another vertical (resp. horizontal) Bell state, see Fig. |6j 
This necessitates the constraint that L is even. We then 
use the fact that for two island Bell states i and j 

X it iX S fl\a,tW,t)j "> XtXjZl^tW^j, (10) 

where A^i acts on the first qubit of Bell state i etc. and 
X\ is a Pauli X on the t qubit of island i. Here Zf 
is Pauli Z on the sign qubit of island or vertex i. The 
island term —AZ^xZi^ — > —AZj, i.e., only acting on 
the t qubit. Note that due to the layout of the sign 
qubits, see Fig. [6j white plaquettes are never frustrated 
since every sign qubit flips the sign of two Ising edges 
of the plaquettes. As every Ising qubit has a sign qubit 
next to it, it implies that a gray square plaquette is frus- 
trated when = Yli^g^Zf = — 1. The parity constraint 
now reads TiiZ\ = 1. The transformed Hamiltonian acts 
with single Zs on the sign qubits, hence for every ba- 
sis state of the sign qubits one obtains a transverse field 
Ising model on the t qubits at the vertices of a square 
lattice on the torus. It is simpler then to represent the 
effect of the sign qubits by associating them with the 
edges of the square lattice. We introduce the gauge bit 
crjj G {—1, 1} for every edge (i,j) which are constrained 
such that H( i j ) Swhitc ^<J(i.j) — 1. One can show that any 
TF Ising model with configuration er satisfying this con- 
straint is unitarily equivalent to a TF model with a sign 
qubit configuration. Thus, modulo unitary transforma- 
tions, we obtain H{a) = A V ( _ rr, r S; .S; - A£\Sf 
where the Ising degrees of freedom are the t qubits and 
the gauge bits are represented by a. The gauge con- 
straint is n^ ^g whitc ^dij = 1 and the parity constraint 
HiSf = 1. Note that the parity constraint makes the 
ground state in the ferromagnetic phase unique, as we 
expect at A = 0. 



We can also study our perturbative expansion, Eq. , 
under this mapping. Each TF Ising model H(a) can 
be handled separately using nondegenerate perturbation 
theory with a finite radius of convergence (^) c . Perturb- 
ing around the (paramagnetic) ground-state | + + . . . +) 
for A = 0, we obtain an effective classical Ising gauge 




FIG. 6: (Color online) The blue horizontal and vertical links 
denote ZZ terms in the Hamiltonian and the black diagonal 
links are proportional to XX on the square-octagon lattice. 
By laying out a basis of Bell states on the blue links, we 
can map the Hamiltonian on the square-octagon lattice to 
a family of signed Ising Hamiltonians in a transverse field 
on a square lattice. The positions of qubit 1 and 2 of each 
Bell state is indicated so that XX links always act between a 
qubit 1 and 2 of different Bell states and we restrict ourselves 
to even L. In the unitarily transformed model, every Ising 
qubit has an additional sign qubit (black dots) next to it. 
If the sign qubit is |1) (marked as red), it flips the adjacent 
red edges from ferromagnetic to anti-ferromagnetic. We can 
alternatively associate a gauge bit (Tjj with every edge 
if we ensure that the product of (Tij around a white plaquette 
is 1. 

theory (modulo some energy shifts), i.e., 
H eB = -A| + + ...+)(+ + .. . + I X 
{& E M n^^e^-J + hi g h er order loops) , (11) 

where | + + . . . +) is the transverse field ground-state. 

Consider what happens in this analysis when the sign 
of the links, Eq. @, is randomly changed (and fixed) 
with respect to the orientation in Fig. [TJ Now we should 
no longer a priori select {c™ = +1} as the ground sector. 
However, all steps can be carried through as long as we 
retain the signs {c™}. A sign change in the links represent 
a sign change in the XX links of the family of Hamiltoni- 
ans on the square-octagon lattice, hence white plaquettes 
can now be frustrated. The proper choice of eigenvalues 
{c^' = ±1} for the ground sector is the one in which no 
square originating from a white plaquette is frustrated. 
This choice may however lead to frustration for the gray 
plaquettes (since C*? -> c™_._ s n (lj)egray ■.). If no 
frustration can be avoided, we may expect the behavior 
of this model to be intrinsically different; this can happen 
on the torus if a single link changes sign. A single link 
sign flip makes — > — for a white and a grey pla- 
quette neighboring the link in the effective Hamiltonian 
H c ff, Eq. (JSJ). But there is no state with lowest possible 
energy — A c ffL 2 . There will be frustration and there are 
10 degenerate toric code ground states with two possible 
excitations all having the same energy of (— L 2 + l)^ 2 -. 
This problem arises since excitations in the toric code 
produce defects on an even number of white plaquettes 
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and an even number of gray plaquettes, hence a single 
sign change of a link does not pick out an excited sec- 
tor as ground sector. For a physical model with open 
boundary conditions such as the surface code, there are 
eigenstates with any pattern of eigenvalues = ±1, 
hence it is possible to avoid frustration. 

Note that in such surface code lay-out, our model will 
have 4 'unlinked' Majorana fermions at the corners of the 
lattice, let us call them c%, C2, C3 and C4. The Hamiltonian 
of this Majorana fermion surface code commutes with 
any string of Majorana operators which starts and ends 
at one of the four corner Majorana fermions. Such strings 
anticommute when they touch the same corner Majorana 
fermion forming a pair of logical operators. Hence the 
spectrum of this Hamiltonian has a double degeneracy 
throughout the entire parameter range; it is the encoding 
of the surface code qubit. When A = 0, the Hamiltonian 
commutes with the four corner Majorana fermions, hence 
the degeneracy is 4 in this regime. This 4-dimensional de- 
generacy is exponentially protected when A is turned on, 
but one of these Majorana qubits (with logical operators, 
say, X = ci, Z — c 2 and Y = icxCz) is only protected by 
fcrmion-parity. At the phase transition, one looses the 
fermion-parity protected Majorana qubit, but keeps the 
Majorana surface code qubit. 

Let us add one last comment on how to handle quasi- 
particle tunneling onto the superconductor island which 
would change the fermion-parity on the island. Such 
event can be viewed as a leakage error with respect to 
our encoding, i.e. an error which maps us outside the 
proper code or ground space. If one can measure the 
parity operator Hq on an island and project the state 
back onto the parity ground space, we reduce such leak- 
age error to no error or a regular error which can either 
be tolerated at low density or be error-corrected actively. 



Projecting Back Lemma 

Let the subspace with fixed eigenvalues {c™} be de- 
noted as H({c™}) and the projector onto this subspace 
as P({c™}). Let H({c™}) be the fermionic Hamiltonian 
after the JW transformation with fixed eigenvalues {c™}, 
i.e the Hamiltonian with links on the square-octagon lat- 
tice, Fig. [Bl We want to consider the spectrum of the 
Hamiltonian H({c™}) on a larger state space only obey- 
ing the parity condition tt = n^C™ = +1 and show that 
the spectrum is identical to that on H({c™}) except for 
additional degeneracies. Let be the projector onto 
the space of states with 7r = UpC™ = +1. We prove 
that all eigenstates of H({c™}) for which P„\ip) ^ have 
the property that P({c?™})\ip) ^ for any {c^} with 
parity n. Since n^C™ commutes with all link operators 
Kj = Xj t \Xj^ (unlike the individual C™), we can con- 
sider the eigenstates in sectors with simultaneously fixed 
eigenvalues of the link operators {Kj} and n. 



Proof: 

Let {Kj} be the set of link operators which commute 
with every H({c%}). Let H({c%})\ip) = E\ip) where 



we can thus assume that Kj 



Consider 



KjP{(c™}M = P{{c'^})K 3 \^) = ±P({c7}M where 
the eigenvalues have changed in sign on two adja- 
cent white plaquettes compared to c™. It is not hard to 
see that there exists a product of link operators K — 
K n . ..K jk such that KP{{c™})\i>) = ±P({c*"})M for 
any pair of {c™} and {c!™} which have the same overall 
parity, i.e., n^c™ = H^c'™ . This is because a link opera- 
tor anticommutes with two plaquette operators C™ and 
hence flips the eigenvalue for two fi, furthermore one can 
generate any pattern of ±1 eigenvalues with fixed parity 
by applying multiple link operators. Since Kj and K are 
unitary, it implies that if some P({c™})\tp) = 0, then all 
P(W™}M = for II^cJ = n M c- and thus P^ip) = 0. 
Hence if P^\ip) ^ 0, then all P({c™})\ip) ^ 0. 
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